Order parameter symmetries for magnetic and superconducting instabilities: 
Bethe-Salpeter analysis of functional renormalization-group solutions 
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The Bethe-Salpeter equation is combined with the temperature-cutoff functional renormalization 
group approach to analyze the order parameter structure for the leading instabilities of the 2D 
t-t' Hubbard model. We find significant deviations from the conventional s-, p-, or d-wave forms, 
which is due to the frustration of antiferromagnetism at small and intermediate t' . With adding 
a direct antiferromagnetic spin-exchange coupling the eigenf unctions in the particle-hole channel 
have extended s-wave form, while in the particle-particle singlet pairing channel a higher angular 
momentum component arises besides the standard d^2_j,2-wave component, which fiattens the 
angular dependence of the gap. For t' closer to t/2 we find a delicate competition of ferromagnetism 
and triplet pairing with a nontrivial pair-wavefunction. 

PACS Numbers: 71.10.Fd; 71.27.H-a;74.25.Dw 



It is by now well established that the superconduct- 
ing order parameter in high-Tc compounds is well de- 
scribed by a {coakx — cos fcy) -momentum dependence: it 
is largest at the Fermi surface (FS) points close to (tt, 0) 
and (0, tt) and vanishes at the FS crossings on the Bril- 
louin zone (BZ) diagonals. Accurate measurements of the 
gap function, however, revealed a slight deviation from 
this d3;2_j,2-wave momentum dependence [1], with a flat- 
ter angular dependence near the nodal points. The sym- 
metry and the details of the momentum dependence of 
the superconducting order parameter are closely related 
to the structure of the effective attractive interaction be- 
tween the electrons. The precise momentum dependence 
of the energy gap therefore contains valuable informa- 
tion about the underlying pairing mechanism, for which 
antiferromagnetic (AFM) spin fluctuations are a viable 
candidate for cuprates [2-5]. 

Another type of unconventional superconductor is the 
layered ruthenate Sr2Ru04 [6], which most likely has 
triplet pairing with p-wave symmetry [7]. It was pro- 
posed that the pairing in this material results from fer- 
romagnetic (FM) spin fluctuations [8,9]. Although in- 
elastic neutron scattering has so far been unsuccessful to 
detect significant low-energy FM spin fluctuations [10], 
this idea finds support from the enhanced tendencies to- 
wards ferromagnetism in the electron doped compound 
Sr2_a;La2;Ru04 [11] and in isoelectronic Ca2Ru04 under 
hydrostatic pressure [12]. 

The important role of AFM and FM fluctuations as 
a possible driving source of singlet and triplet supercon- 
ductivity, respectively, was emphasized early on in the 
pioneering work in Refs. [13-15]. Recently, the inter- 
play of antiferromagnetism and d-wave superconductiv- 
ity (dSC) and ferromagnetism and p-wave superconduc- 
tivity (pSC), respectively, was reconsidered within the t- 



t' Hubbard model using functional renormalization-group 
(fRG) techniques [16-20]. Early versions of fRG [16-18], 
which used the momentum cutoff procedure, were unable 
to search for ferromagnetism. This drawback is overcome 
in the temperature-cutoff fRG approach (TCRG) [19], 
which proved successful in describing both, AFM and FM 
instabilities together with singlet- and triplet supercon- 
ducting pairing in the weak-coupling t-t' Hubbard model 
and its extensions in a broad parameter range [19-21]. In 
the previous fRG analyses it has been a common prac- 
tice to assume order parameter structures, which have 
the form of square lattice basis functions with s-, p- or 
d-wave symmetry [16-21]. Although it was recognized 
that in the RG flow of the order parameter susceptibili- 
ties also the momentum dependence acquires specific cor- 
rections to their initial form [17,19-22], these corrections 
were so far not analyzed in detail. Indeed, actual order 
parameter structures may have admixtures of different 
symmetry components, and it is necessary to develop an 
unbiased method to determine their precise momentum 
dependence. 

In the present paper we use the Bethe-Salpeter equa- 
tions to extract eigenf unctions and eigenvalues of the 
effective interaction in the particle-particle (pp) or the 
particle-hole (ph) channel - similarly as in previous quan- 
tum Monte Carlo studies [23]. Here we consider a com- 
bination of the Bethe-Salpeter equations and the fRG 
approach. We choose the TCRG version [19] as the most 
suitable tool, because the effective model obtained within 
this scheme does not contain any unintegrated degrees of 
freedom even at the intermediate stages of the RG flow. 

We apply this procedure to identify the order param- 
eter structure for the leading instabilities of the 2D t-t' 
extended Hubbard model H^^ = H — (p — 4:t')N with 
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where Uj = t for nearest neighbor (nn) sites i and j and 
= —t' for next-nn sites (i, t' > 0) on a square lattice; 
for convenience we have shifted the chemical potential 
/i by At' . In Eq. (1) we have included a direct nn spin 
exchange interaction J; Si = c[^craf}Cii3/2, and a de- 
notes the Pauli matrices. While such an interaction is 
generated from the on-site Coulomb repulsion at strong- 
coupling, we add it here as an independent interaction 
in the weak-coupling regime, where the RG scheme is 
applicable. 

We follow the many-patch fRG version for onc-particle 
irreducible Green functions as proposed in Rcf. [19]. This 
TCRG scheme uses the temperature as a natural cutoff 
parameter and accounts for excitations with momenta 
far from and close to the FS, which is necessary for the 
description of instabilities arising from zero-momentum 
ph scattering, e.g. ferromagnetism. Neglecting the fre- 
quency dependence of the vertices, which is expected to 
have minor relevance in the weak-coupling regime, the 
RG differential equation for the interaction vertex has 
the form [19] 
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dT 



dT 
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where o is a short notation for summations over interme- 
diate momenta and spin. 
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and fris) is the Fermi function. The upper sign in 
Eq. (3) is for Lp/j and the lower sign for Lpp, respec- 
tively. Eq. (2) has to be solved with the initial condition 
Vjo (ki, k2, ka, k4) = C/; the initial temperature is chosen 
as large as Tq = AQQt. 

We discretize the momentum space in Np ~ 48 patches 
using the same patching scheme as in Ref. [19]. This re- 
duces the integro-differential equations (2) and (4) to a 
set of 5824 differential equations, which were solved nu- 
merically. The evolution of the vertices with decreasing 
temperature determines the temperature dependence of 
the susceptibilities according to [19,24] 
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FIG. 1. Diagrammatic representation of the three-point 
vertex TZ^'^ (solid triangle), which has the structure of a ver- 
tex diagram in a static external field acting in the channel m 
with the r-th initial vertex function /^'^' defined in Eq. (6) 
and shown as an open triangle. The specific example shows 
the phsQ channel. The circle represents the reducible elec- 
tron-electron interaction vertex. 



m denotes one of the possible channels: ph spin (phs), 
which traces FM and AFM instabilities, ph charge (phc) 
for the Pomeranchuk instability (PI) [25] or d-density 
wave (DDW) [26] , or pp for superconducting singlet and 
triplet instabihties; q„i = Q =(7r,7r) for AFM and DDW 
and qm = otherwise. In the following we use the no- 
tation rnqm to denote the specific channel. The index r 
in Eq. (4) represents the symmetry of the corresponding 
channel. The upper signs and pp-indices in Eq. (4) refer 
to the superconducting instabilities (dSC and pSC), and 
the lower signs and ph-indices refer to the charge and 
magnetic instabilities. 

The three-point vertices TV^'^ can be considered as 
vertices describing the propagation of the electron in a 
static external field; their diagrammatic representation is 
shown in Fig. 1. The initial conditions at Tq for Eqs. (4) 

(r) 



are 



fk^ and Xr^ 



0, where the function /j^' 



belongs to one of the irreducible representations of the 
point group of the square lattice {Aig or -Bi.2g), e.g. 



Ar) 



A- 




cos ky d. 



y2-wave (Big), 
p-wave {B2g), (6) 

s-wave (^ig) 



with a normalization coefficient A = (l/A^)X]k/k- 
momentum dependence of the vertices TZ^'^ changes dur- 
ing the RG fiow and at low temperatures it is expected to 
refiect the momentum dependence of the ground-state or- 
der parameter of the selected symmetry. However, these 
vertex functions depend on the initial choice of the func- 

(r) 

tions , and therefore they can not serve as an un- 
biased tool to obtain the structure of the ground-state 
order parameters. 

To perform an alternative analysis, which is not based 
on the knowledge of the wave functions at high temper- 
atures, we consider the solution of the Bethe-Salpeter 
equations [23] 



^rpft(k,p;p,k)Lph(p,p-f q™)0P'' = \ph(t>\ 



ph 
k ' 



yT(k,k',k' + q„) phs (AFM, FM), 
yT(k,-k + q„„k') pp(dSC, pSC). 



- J2 rpp(k, p; p, k)ipp(p, p)rj' = AppC , (7) 



where Tph and Tpp denote irreducible vertices in ph and 
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pp channels, respectively. Exploiting the connection to 
reducible vertices, Eqs. (7) can be rewritten as [23] 



A, 



1-A 



m.l 
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where I enumerates eigenvalues and -functions for a given 
channel m. The reducible vertices (k, p) can be di- 
rectly extracted from the fRG flow according to Eq. (5). 

The value Xm,i = 1 corresponds to an ordering in- 
stability with the symmetry of the eigenfunction 0™' . 
Therefore, tracing the temperature dependence of eigen- 
values and -functions allows to identify both, the lead- 
ing instabilities and their concomitant order parame- 
ter structure. We stop the fRG flow at the tempera- 
ture Tx, where the maximum interaction vertex Knax = 
max{F(ki,k2;k3,k4)} reaches the value 20t. Tx should 
be understood as a crossover temperature into a renor- 
malized classical regime with exponentially large correla- 
tion length [21], we have verified that the following results 
for the eigenfunctions are only weakly dependent on the 
choice of Tx ■ 

Eigenfunctions and eigenvalues of the Bcthc-Salpctcr 
equations provide more detailed information than can 
be obtained from the momentum dependence of the 

three-point vertices TZ^'^, which determine the order- 
parameter susceptibilities according to Eq. (4). Unlike 
the vertices TZ^'^, the solutions of the Bethc-Salpeter 
equations do not depend on the initial choice of the func- 

(r) 

tions and, therefore, do not require a priori the 
knowledge of the symmetry of the leading instability. 
The three-point vertices 



/r+Er^(k,p)iJ/.,pp(p,±P + q™)/^' (9) 



described by the diagram of Fig. 1 can be obtained di- 
rectly from the solutions of the Bethe-Salpeter equations. 

(r) 

Expanding the functions in terms of eigenfunctions 
of the Bethe-Salpeter equation and using Eqs. (8), we 
obtain 



i'm,l 



(10) 



The vertices 7^™''^ found in this way, however, do not 
necessarily coincide with those obtained directly from the 
RG procedure, Eq. (4), because generally the one- loop 
approximation does not exactly reproduce the solution 
of the corresponding Bethe-Salpeter equation. The only 
case when the three-point vertices TZ^''^ and TV^'^ co- 
incide is the ladder approximation when either only the 
pp-channcl or one of the ph-channels in Eq. (2) is re- 
tained. Nevertheless, as we will see below, the results for 
the vertex TZ^'^ from Eqs. (4) are qualitatively similar 
to those found from Eq. (10). 

It is clearly established from Eq. (10) that TZ^'"^ 
(and similarly Ti^'^) is actually a mixture of different 
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FIG. 2. Eigenvalues (a) and angular ^kj^ -dependence on 
the FS of the eigenfunctions 0^'' = 0; of the Bethe-Salpeter 
equation and the normalized three-point vertices TZ^'^ (de- 
noted as n-m,r) at T = Tx (b,c,d) for t' = O.lt, U = 2t, 
J = H = 0. Tx is the lowest temperature reached in the fRG 
flow (see (a)). 



eigenfunctions (f)^' , whose weights arc proportional to 
1/(1 — Xm,i) with coefficients determined by the overlap 

(r) 

of the eigenfunction with /^/. If one of the eigenval- 
ues X„, i is much closer to unity than all the others, the 
corresponding term in the sum over I in Eq. (10) is ex- 
pected to dominate and TZ^'^ essentially coincides with 
the corresponding eigenfunction ^J^''. 

Below we discuss the results of the ntunerical solution 
of the Bethe-Salpeter equations and compare them to the 
results for the three-point vertices TZ^'^ for different pa- 
rameter regimes. We start in Fig. 2 with results at the 
van Hove (vH) band filling (/i = 0) for small t' = O.lt, 
J = 0, and U = 2t. The AFM phsQ instability has the 
largest eigenvalue at Tx in agreement with previous fRG 
work based on the analysis of order parameter suscepti- 
bilities [19]. The corresponding eigenfunction (Fig. 2c) 
has a slight variation around the FS with an enhancement 
near (tt, 0) and (0,7r), which most likely originates from 
the vH singularity nature of these points. The eigenfunc- 
tions in the subleading phcO and ppO channels are sizable 
near (tt, 0) and (0, tt,) only. Although the eigenvalue of 
the zero-momentum ph instability in the charge channel 
(phcO) is relatively small at T = Tx (A ~ 0.5), it rapidly 
increases at low temperatures. The corresponding eigen- 
function (Fig. 2d) is antisymmetric with respect to 90° 
rotation, thus providing the possibility for the Pomer- 
anchuk instability with a spontaneous d-wave like defor- 
mation of the FS [27,28]. The eigenvalues for the PI at 
t' = O.lt are larger than those at = (not shown), 
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FIG. 3. Same as in Fig. 1 for = 1.5t, J = 0.3t 



FIG. 4. Same as in Fig. 1 for = O.St, U = 2t, J = 0.3t. 



which supports the conclusion of Rcf . [28] , that a finite 
value of t' at vH band fillings enhances the tendency to- 
wards a Pomeranchuk instability. 

The momentum dependence of the three-point vertices 
is qualitatively similar to the leading eigenfunctions 
in the ppO and phsQ channels, however with a smaller 
variation around the FS. In the phcO channel we observe 
a stronger difference of the vertex T^p'^^.d-wave^ which has 
almost the da,2_y2-wave form, from the corresponding 
eigenfunction of the Bethe-Salpeter equation. 

To investigate the role of t'-induced frustration of an- 
tiferromagnetism, we turn on a small direct exchange in- 
teraction J = 0.3t to strengthen the AFM correlations, 
and decrease simultaneously the value of f/ to 1.5t in or- 
der to remain in the weak-coupling regime. Fig. 3 shows 
the resulting changes. At finite J the eigenfunctions in 
the ppO and phcO channels are essentially nonzero all 
around the FS, however they arc flatter than the dx2_y2- 
wave function. The eigenfunction in the phsQ chan- 
nel is of extended s-wave form. The eigenvalue for the 
phcO channel is smaller than at J = 0, implying that 
stronger antiferromagnetism weakens the tendency to- 
wards a Pomeranchuk instability. Similarly to the J — 
case, we observe a smaller variation of the vertices TZ^'^ 
around the FS, than in the momentum dependence of 
the corresponding eigenfunctions of the Bethe-Salpeter 
equation. 

With increasing t' to OM (Fig. 4) the largest eigen- 
value occurs in the singlet dSC channel. The ppO pair 
wavefunction maintains its shape with a slight deviation 
from the (cos — cos fcy)-form, i.e. a flattening near the 
BZ diagonal. In the phsQ and phcO channels we observe 
again a weaker momentum variation of TiJ^'^ . The ver- 



tex ^^j^ instead is very close to the shape of the 

corresponding Bethe-Salpeter eigenfunctions in the ppO 
channel. 

For t' = 0.45t a ferromagnetic instability is expected 
[19,20], and we start again with the vH band fllling case 
for /i = 0. We increase the interaction strength to f7 = 3t, 
since the corresponding crossover temperatures for the 
FM instability are lower. The largest Bethe-Salpeter 
eigenvalues arise in the FM (phsO), pSC (ppO triplet) 
and AFM channels (phsQ). Remarkably, the wavefunc- 
tion in the triplet ppO channel deviates signiflcantly from 
the sin form (see Fig. 5b) , but its eigenvalue remains 
smaller than the eigenvalue for ferromagnetism. This de- 
viation is also well reproduced by the corresponding mo- 
mentum dependence of the vertex t^pp^'P"™™^ 

On moving away from the vH band filling at t' = 0A5t, 
we observe a further increase of the eigenvalue in the 
triplet pairing channel (Fig. 6), which however remains 
smaller than the eigenvalue of the ferromagnetic insta- 
bility at the lowest temperature we can safely reach in 
the fRG flow. Simultaneously, the eigenfunction in the 
triplet superconducting channel slightly distorts towards 
the wave function y^P"'^'*™) ^ strong deviations persist 
in both, the eigenfimction of the Bethe-Salpeter equa- 
tion and the vertex -^ppO'P-™'*™ Singlet superconductiv- 
ity, which was not among the dominant instabilities at 
/U = 0, is also significantly enhanced for /j, > 0. There 
are two eigenfunctions, which are symmetric and anti- 
symmetric with respect to reflection at the BZ diago- 
nal with almost equal eigenvalues in the pp singlet chan- 
nel. Both eigenfunctions are essentially nonzero at (tt, 0) 
and (0,7r) only, and therefore not expected to be stabi- 
lized thermodynamically. Because of the presence of two 
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nearly-degenerate eigenf unctions, the three-point vertex 
y^ppo.d-wave ^j^-^ substantially deviates from the 
momentum dependence of the leading eigenfunction. 

In conclusion, we have investigated the symmetry of 
the leading instabilities of the 2D t-t' Hubbard model us- 
ing as a novel tool the combination of the Bethe-Salpeter 
equation and the fRG approach. Although in most cases 
the leading instabilities coincide with those from suscep- 
tibility based analyses, the true shape of the eigenfunc- 
tions significantly differs from the s-, p- or d-wave basis 
functions. At small t' we find in particular a flattening 
of the eigenfunction in the pairing ppO channel near the 
nodes - in qualitative agreement with experimental data 
for cuprates [1]. The addition of a direct spin exchange 
interaction to the Hubbard model at weak-coupling was 
essential to reproduce a (cos fc^ — cos A;y)-like form of the 
superconducting gap at low and intermediate t'. The in- 
stability towards triplet pairing at larger t' also shows 
a substantial deviation from the standard p-wave sin kx~ 
form. In many cases we have investigated, the eigen- 
functions of the Bethc-Salpcter equations strongly de- 
viate from the corresponding three-point vertices TZ^'^ , 
which enter the order parameter susceptibilities, due to 
either a qualitative difference of the cigenfunctions from 
the basis functions j^^^'P''^)-'"^'^'^ or the presence of sev- 
eral almost degenerate cigenfunctions. The proposed new 
technique may prove most useful in future studies of mag- 
netic and superconducting instabilities without presup- 
posing a special momentum dependence of the candidate 
order parameters. 
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